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ABSTRACT 



I 1 I — 

Following iFord et all (|2011 



20121) and 



1 I 1 

Steffen et al.l (|2012aj) we derived transit 



timing variations (TTVs) of the Kepler KOIs, using the first twelve quarters of 
the Kepler data. Simultaneously with its timing, we obtained for each transit 
its duration and depth, presenting the results as a catalog for the community 
to use. We derived a few statistics of our TTVs that could be used to indicate 
significant variations. Including systems found by previous works, we pointed 
out 100 KOIs that showed highly significant TTV, 9 of which had short-period 
TTV modulations with small amplitudes. We considered two effects that could 
cause apparent periodic TTV — the finite sampling of the observations and the 
interference with the stellar activity, stellar spots in particular. We briefly discuss 
some statistical aspects of our detected TTVs. We show that the TTV period is 
correlated with the orbital period of the planet and the TTV amplitude. 



Subject headings: planetary systems - planets and satellites: detection - techniques: 
miscellaneous 
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Introduction 



Since 2009 May 2, the Kepler spacecraft has been collecting scienc e-quahty photometric 



data of more than 150,000 stars. Based on the first 5 months of data, 



Borucki et al 



(12011 



hereafter Bll) identified 1235 planet candidates as sociated with 997 host stars. Analysis 



of the first 16 months of data ( Batalha et al. 



20121 . hereafter B12) yielded an additional 



1091 viable planet candidates — termed as Kepler objects of interest, or KOIs. The 
almost uninterrupted accurate Kepler light curves of these KOIs enable the community to 
detect minute changes in the observed transit light curves. This is especially true for the 
individual times of transit, which for some KOIs show variation (=TTV) relative to a linear 
ephemeris that assumes a constant Keplerian orbit. These TTVs can indicate a dynamical 
interaction with addition al obj ects in the s ystem , as was predicted by the seminal works of 



Holman fc Murray 



fl2005[ l and 



Agol et al. 



( I2OO5I ). Indeed, TTVs turn out to be a crucial 



tool i n the study of systerns with known multip 



2010 



Lissauer et al. 



2011a : 



Cochran et al. 



2011 



e transiting planets ( e.g.. 



Fabrycky et al.li201 



Holman et al. 



Steffen et al. 



2012bl ) 



However, TTVs can do much more. They may indicate dynami cal interactions wit 
unseen, otherwise undetected, additional objects in the system (e.g 



Ballard et al 



20111) 



Therefore, i t can be us e: 
the work of Ford et al. 



hereafter F12) and 



ul to perform a systematic search of all KOI s, as was d o ne in 



2011 



he reafter Fll) and was continued with 



Ford et al. 



Steffen et al. 



(2012 



(l2012al . hereafter SI 2) work based on the first six quarters 
of Kepler data. This paper is a follow-up of Fll, F12 and S12 studies (see also the catalog 
of Rowe et al., private communication), presenting a systematic analysis of the first twelve 
quarters of the Kepler data of all KOIs. The goal is to produce an easy-to-use catalog that 
can stimulate further analysis of interesting systems and statistical analysis of the Kepler 
KOIs with significant TTV sample. 



After presenting the details of our pipehne and the catalog itself (Section [2]), we derive 
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a few statistical characteristics of each TTV series that can identify the ones with significant 
variations (Section 3). Sections 4 & 5 hst 100 systems with highly significant TTVs, and 
Section [H] comments on some interesting systems, in particular the ones for which the 
derived TTVs could be of non- physical origin. In Section [7] we present a few basic statistical 
features of the sample of the 100 systems, and briefly discuss the possible use of the catalog. 



2. Analysis of the transit light curves 

The Kepler long-cadence data were downloaded from the archive.stsci.edu ftp server 



(ftp://archive.stsci.edu/pub/kepler/lightcurves/tarfiles). We obtained data for 2318 KOIs 



listed in Bll and/or B12 ( http : / /archive. stsci.edu/kepler/ planet-p andidates.html). We 
did not analyze 188 KOIs: 

• 51 KOIs of SNR < 10, which we opted not to analyze. 

• 21 KOIs did not have enough information for our pipeline. 20 KOIs had only one 
transit and one KOI did not have a measured transit duration. 

• 116 KOIs had very short transit duration of less than 1.5 hours, which resulted in less 
than 3 points per transit. 

Out of the remaining 2130 KOIs, our program did not succeed in analyzing 7 KOIs: 
5 KOIs produced a bad model and 2 had very long orbital periods resulting in only 2-3 
measurements, which left us with 2123 KOIs that are included in Tables [1] & [2l However, 
86 KOIs had less than 7 acceptable TTV measurements, and therefore we did not include 
them in the statistical analysis. We were left with 2037 KOIs (see Table [3]). 

We started by folding the Kepler light curve of each KOI with the ephemeris of Bll 
and/or B12, in order to obtain the best possible template for the transit light curve (see 
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Fll for de tails). In order t o res erve some flexibility for the fitted model we opted not 



to use the iMandel &: Agoll (|2002| ) model, which used a specific transit shape with a given 
limb-darkening function, ignoring, for example, any stellar gravity darkening. We used 
instead an analytical function that could mimic different transit shapes. In particular, w e 



allow ed a slig ht asymmetry in ou r model, based on the example of KOI-13 ( ISzabo et al 



2011 



see also 



Mazeh et al. 



2OI2I ). Such a flexibility came with a price. When the folded 



light curve of a transit had a low SNR, the algorithm adopted a strange shape which was 
obviously wrong. For such cases we added an option to fit the data with a second model, 
which is less flexible but maintained a realistic transit profile. 

We used the best-fit transit model as a template to measure the actual timing of each 
individual transit (=TT) and derive its TTV — the difference between the TT and the 
expected time, based on the linear ephemeris. The TT derivation was performed while 
allowing the duration — TDV, and depth — TPV, of each transit to vary as well. As our 
template was just a mathematical function, the TTV, TDV and TPV were equivalent to 
moving the center of the template or stretching it in the time and flux dimensions. In our 
approach, we searched for the minimum of the sum-of-squared residuals, the standard 
function, in the three parameter space. 

Although the main emphasis of this paper was the TTVs of the KOIs, we opted to vary 
the three parameters of the template simultaneously, as we found a few KOIs for which the 
transit duration or depth did vary significantly, either because of physical proc esses, or as a 
results of some observational effects. One example is, again, KOI-13, for which 



Szabo et al. 



(l2012aj) have found some indications for a long-term variation of the impact parameter. 



which is equivalent to the variation of ^ t he tran sit duration. Our analysis, based now on 
twelve quarters, confirmed iSzabo et al.i (i2012a) results and are presented in Figure [H One 
can see in that figure the highly significant linear duration variation of KOI-13, which 
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amounts to ~ 1% peak-to-peak modulation over the entire span of the data. For such cases 
the simultaneous analysis of TTV, TDV, TPV was an advantage, and, in principle, could 
yield better TTVs. 



KOI 13.01 TDV 



-0.005 





Median Errors 












'% ' * ' , 






-'V- C •• '/ 











Time (BJD -2454900) 

Fig. 1. — The TDV (duration variation) of KOI-13.01. Each point represents our best 
estimate for the deviation of the transit duration from its averaged value. 

However, for low-SNR transits the three-parameter simultaneous minimization of the 
function could give completely erroneous results, based on some local minimum in a 
noisy surface. We therefore had to strictly vet our results and reject any suspicious 
findings. As a result we did not report in our catalog TTV, TDV or TPV for some transits. 
We believe that many of those transits would yield suspicious TTVs anyway, and therefore 
the price we paid for the simultaneous fitting was not too high. 



2.1. The model for a transit 

We select our best model for the folded transit among two templates — a Legendre- 
based and a Fermi-based model. Both models assume that the light curve has a unity 
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value outside the transit, which lasts from ti until ^4. The choice between the two is made 
according to their goodness of fit, i.e. the model with the lower x^- We further make sure 
that the model has one local minimum only. 

The Legendre-based model has the form: 

N 

F{t) = J2 i^kL'' (r) + SkSin {nk/2 ■ r) + Ckcos {7ik/2 ■ r)] , (1) 

k=l 

where L'' is the Legendre polynomial of order k, r is the normalized phase of the transit, 
such that at the beginning of the ingress t {ti) — —1, and at the end of the egress 
T {ti) — +1, ^fe, Sk and Cfe are linear parameters found through a linear optimization model, 
and A'^ = 10 is the maximum order we allowed. We optimize the model through varying 
the phases of ti and within the orbital period. We avoid local bumps in the model by 
reducing its order and by using local linear fits when needed. 

The Fermi-based model has the form: 



F (r) = 1 + M 



g(r+<^+M)/s _|_ 1 ~'~ g(T+¥'-;i)/s _|_ ^ 



(2) 



where r is the phase of the transit, as for the Legendre model, and 99, s, n and M are 
free parameters which define the phase shift, steepness, width and depth of the transit. 
In order to obtain a more "transit-like" shape, the points between the ingress and egress 
are replaced by a parabola, with its width as another free parameter, under the analytical 
constraint that the resulting function and its first derivative are both continuous. We used 
Markov chain Monte Carlo (MCMC) process to find the best values of the five parameters. 

In most KOls the selected function was the Legendre-based model, while the 
Fermi-based one was typically selected only when the SNR of the folded light curve was 
low. 
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2.2. Finding the best TTV, TDV and TPV 

The derivation of the TTV, TDV and TPV was done after we fitted a polynomial 
to the light curve on the two sides of each transit, in order to remove the stellar and 
instrumental long-term photometric variations at the time of each transit. We minimized 
the standard function, assuming each cadence has the same error. The errors on each 
TTV, TDV and TPV were estimated by searching for the first point for which A^^ = 1 
relative to its minimum value. The scatter around the polynomial fit before and after each 
transit was used as our best estimate for the error of the individual Kepler measurements 
at that transit. 

For each KOI, we rejected outlier TTVs, either because the TTV value was too 
different from the other TTVs of this KOI, or because the TTV error estimate was too large. 
Usually a large err or meant that som e photometric measurements during that transit were 



erroneous. Like in 



Ford et al. 



(120111 ). we iterated the procedure, aligning the transits based 
on the measured TTVs, in order to generate a better transit model, and then recomputed 
the TTV, TDP and TPV. 

Table [T] lists the modified ephemerides of the KOIs, based on our analysis, and the 
duration and depth of their transit, derived from the folded light curve. The transit duration 
is given as a fraction of the orbital period and the depth is given in units of the stellar 
intensity outside the transit. Table |2] lists our derived TTV, TDV and TPVs for 170438 
transits. We list the transit number, its expected time, and then its TTV/TDV/TPV with 
their uncertainties, respectively. The TDV and the TPV are given in units of the transit 
duration and depth. 
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Table 1: Linear ephemerides of the KOI transit s, together with their durations and depths 



KOI 




To- 
[d] 




Period 
[d] 


Duration*^ 


Depth^ 


1.01 


55, 


.762503 


2 


.47061339 


0. 


.0369 


0.01412 




±0, 


.000008 


±0, 


.00000003 








2.01 


54, 


.357917 


2 


.20473527 


0. 


.0837 


0.00667 




±0, 


.000018 


±0, 


.00000006 








3.01 


57, 


.812653 


4, 


.88780170 


0. 


.0242 


0.00433 




±0, 


.000033 


±0, 


.00000027 








4.01 


90, 


.526022 


3 


.84937091 


0. 


.0334 


0.00131 




±0, 


.000311 


±0, 


.00000188 








5.01 


65, 


.973225 


4, 


.78032754 


0, 


,0215 


0.00096 




±0, 


.000232 


±0, 


.00000179 








7.01 


56, 


.611798 


3 


.21366637 


0, 


,0615 


0.00073 




-1-0 


.uuuooy 




0000090"^ 
.uuuuuzuo 








10.01 


54, 


.118496 


3 


,52249844 


0, 


.0425 


0.00913 




±0, 


.000055 


±0, 


.00000031 








12.01 


79, 


.596487 


17, 


.85521360 


0. 


.0208 


0.00901 




±0, 


.000344 


±0, 


.00000992 








13.01 


53, 


.565129 


1 


.76358761 


0. 


.0858 


0.00457 




±0, 


.000009 


±0, 


.00000003 








17.01 


54, 


.485703 


3 


.23469926 


0. 


.0514 


0.01057 




±0, 


.000032 


±0, 


.00000018 









"BJD-2454900. 

''In units of the orbital period 

^In units of the stellar intensity outside the transit 
(This table will be available in its entirety in a machine- readable form. A portion is shown here for guidance 
regarding its form and content.) 
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Table 2: TTV, TDV and TPV of the Kepler KOIs 



n 




TTV„ 




TDV„ 


CTn 


TPV„ 


CTn 




[d] 


[min] 


[min] 













55.7625 


-0.034 


0.096 


0.0055 


0.0027 


-0.0116 


0.0026 


1 


58.2331 


0.06 


0.077 


0.0025 


0.0021 


-0.0129 


0.0021 


2 


60.7037 


-0.13 


0.1 


0.0045 


0.0026 


-0.0149 


0.0028 


3 


63.1743 


-0.03 


0.13 


0.0018 


0.0032 


-0.016 


0.0034 


4 


65.645 


-0.21 


0.11 


0.0015 


0.0026 


-0.0032 


0.0029 


5 


68.1156 


-0.021 


0.087 


-0.0007 


0.0022 


-0.0046 


0.0024 


6 


70.5862 


0.06 


0.1 


-0.0016 


0.0027 


-0.0088 


0.0028 


7 


73.0568 


0.13 


0.11 


0.018 


0.0032 


-0.0333 


0.003 


8 


75.5274 


0.17 


0.11 


0.0043 


0.0034 


-0.0045 


0.0032 


9 


77.998 


0.04 


0.11 


-0.0056 


0.0035 


0.0027 


0.0032 



"BJD-2454900. 

(This table will be available in its entirety in a machine-readable form. A portion of the table of KOI- 1.01 
is shown here for guidance regarding its form and content.) 
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3. Identification of KOIs with significant TTV 

As the main interest of this study is the TTVs of the KOIs, the next sections 
concentrate on the analysis of the derived TTV. Analyses of TDVs and TPVs are deferred 
to a later paper. 

In order to identify KOIs with significant TTVs, we compute a few statistics (see 
Fll, F12 and S12) that characterize some modes of the scatter of the derived TTVs, all 
presented in Table [3l 

• For each KOI, we compare the MADtj-v — the median absolute deviation of the 
TTVs, which is a measure of their scatter, and (Tttv, their median error (see Fll 
for details). High values of MADttv, relative to <Jttv, may indicate significant 
TTV, especially because the MAD statistic is not so sensitive to outliers. However, 
the derived ratio relies on our estimate of the TTV error, which by itself depends 
on the estimated error and the nature of the noise of the Kepler measurements. As 
we are not sure what contribution red noise makes to the total noise of the Kepler 
measurements, we do not rely on our estimate of the TTV error, and therefore do not 
use the MAD to a ratio to identify KOIs with significant TTVs. The two values are 
given in the table for completeness. 

Another drawback of the MAD to cr ratio is its insensitivity to the order of the 
residuals. That is, any permuta tion of the residuals yields the same two v alues. However, 



as pointed out by Fll (see also lAgol et al. 



2005 



Holman &: Murray 



20051 ). the expected 



timescale of the dynamical interaction between planets is in most cases larger than the 
orbital period of the transiting planet, and therefore we expect long-term correlation in the 
planet's TTVs, if indeed the planet is subject to a dynamical perturbation. 

We therefore add three statistics that can indicate long-term correlation of the TTVs. 
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The Lomb-Scargle (LS) periodogram (e.g., S12), which searches for a cosine-shape 
periodicity in the series of TTVs. We identify the highest peak in the periodogram and 
assign a false-alarm probability (FAP) to the existence of the associated periodicity 
in the data. This is done by calculating similar 10^ LS periodograms with different 
random permutations of the same TTVs, and obtaining the highest peak in each of 
these periodograms. The best period and log of its p- value are listed in Table El 

A long-term polynomial fit to the series of TTVs. A good polynomial fit usually 
indicates a long-term modulation with a time scale longer than the data span. We 
search for a polynomial with degree lower than 4, choose the best fit and check its 
significance with the T test (e.g., Fll). The log of the p- value of the best polynomial 
fit is also given in Table O 



The alarm score A of 



Tamuz. Mazeh &: North 



(120061 ) . which is sensitive to correlation 



between adjacent TTVs. The value of A refiects the number of consecutive TTVs 
with the same sign, with o ut as suming any functional shape of the modulation (see 



Tamuz. Mazeh fc North 



(120061 ) for a detailed discussion). This is done relative 
to a constant, zero TTV assumption. We assign a false- alarm probability to the 
occurrence of the obtained score by calculating alarm scores for 10^ different random 
permutations of the same TTVs. The alarm score and the log of its p-value are listed 
in Table El 



Table El can be used to identify KOIs with significant TTVs of various timescales. 



4. KOIs with significant TTVs 



In this section we point out 91 systems with highly significant TTVs, plotted in 
Figures [214T31 and summarized in Table HI Seventy six of those KOIs show some periodicity. 
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Table 3: Statistical parameters for TTVs of Kepler KOIs 



KOI 


(T 


MAD 


LS 


LS 


p-LS 


A 


p-^ 


Pol. 










Period 


Peak 








Deg. 






[min] 


[min] 


[d] 




[log] 




[log] 




[log] 


1.01 


0.0853 


0.0628 


10.8 


5.43 


-0.1 


-0.033 


-0.5 


2 


-0.3 


2.01 


0.2392 


0.1556 


5.36 


6.05 


-0.1 


0.233 


-1.2 


2 


-0.3 


3.01 


0.1735 


0.1802 


34.8 


5.42 


-0.2 


-0.499 


-0.1 


2 


-0.3 


4.01 


2.2182 


2.031 


10.41 


6.2 


-0.4 


-0.419 


-0.1 


2 


-0.3 


5.01 


1.6114 


1.0833 


14.64 


6.62 


-0.6 


-0.183 


-0.3 


1 


-0.3 


7.01 


3.3156 


2.2459 


45.48 


5.59 


-0.1 


-0.066 


-0.5 


1 


-0.3 


10.01 


0.5838 


0.4043 


75.82 


5.36 


-0.1 


-0.2 


-0.1 


2 


-0.3 


12.01 


0.6881 


0.9077 


1061.48 


7.64 


-2 


0.873 


-1.6 


2 


-0.4 


13.01 


0.1689 


0.0921 


5.72 


6.18 


-0.1 


-0.244 


-0.1 


3 


-0.3 


17.01 


0.3259 


0.2713 


10.7 


9.04 


-1.4 


-0.102 


-0.3 


2 


-0.3 



(This table will be available in its entirety in a machine-readable form. A portion is shown here for guidance 
regarding its form and content.) 
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with time scales ranging from ~ 200 to ~ 1400 days and amplitudes in range of ~ 10 to 
~ 1200 min. For each of the 76 systems, we suggest a possible continuous fit to the TTVs, 
composed of a straight line, which can present a correction to the orbital period of the 
transiting planet, together with a cosine function with the best-found period. For 47 KOIs, 
the period, its error and amplitude are given in the table. For 28 systems the period found 
was too long and we can not derive its uncertainty. In those cases, the period given is 
just an approximation. In one special case — KOI-142.01, we fit a straight line with two 
different cosine functions. 

For 11 KOIs, the TTV series do not include a maximum and a minimum, and therefore 
a cosine function can not be fitted, probably because the time scale of the modulation is 
longer than the time span of the data. In those cases, we fit the TTV with a long-term 
parabola only, and a note is written in the plot and in Table IH For four systems — 
KOI-1241.01, 1285.01, 1452.01 and 1546.01, neither a cosine function nor a parabola could 
be fitted, but the TTVs look nevertheless significant (see Section [6] for s short discussion). 

Table H] lists the KOI number, the orbital period of the transiting planet and the 
model we use, either a Cosine function, "C", or a polynomial "P". For a Cosine fit, the 
TTV period and its error, when available, and amplitude are given too. Also listed are 
the number of TTV measurements, the number of KOIs in this system and references to 
previous studies, when available. The table includes references to Section [6|, where we 
briefly comment on some of these systems. 
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Table 4: KOIs with significant TTV 



KOI 


Period 


Model 


Period 




Amp 




N 


Multi- 


Ref. 










[dl 


[mill] 


[min] 




plicity 




42.01 


17.83 


c 


990 




13.95 


0.86 


51 


1 




84.01 


9.29 


C 


309 


31 


4.47 


0.41 


101 


1 


Kcplcrl9b 


103.01 


14.91 


c 


261.7 


8.1 


25.64 


0.83 


62 


1 




137.01 


7.64 


c 


262 


22 


4.31 


0.41 


92 


3 


1''* Kcplci-18c 


137.02 


14.86 


c 


266 


25 


3.93 


0.38 


53 


3 


i'2'4Keplerl8d 


*142.01 


10.95 


c 


632 


38 


669 


15 


88 


1 


1,2 








349 


19 


109 


3 








148.01 


4.78 


c 


404 


68 


4.18 


1 


144 


3 


^ Kepler48b 


148.02 


9.67 


c 


411 


64 


3.46 


0.8 


74 


3 


^ KeDler48c 


152.02 


27.40 


c 


640 


320 


17.6 


3.3 


24 


3 


13 


157.04 


46.69 


c 


590 


130 


24.9 


4.7 


18 


6 


^'^ Kepler 1 If 


168.01 


10.74 


c 


491 


78 


22.2 


3.6 


69 


3 


^'■^ Kepler23c 


*190.01 


12.26 


c 


271 


18 


4.14 


0.33 


65 


1 


16 


227.01 


17.70 


c 


1600 




1042.3 


6.9 


44 


1 


1>2 


244.01 


12.72 


c 


362 


45 


1.04 


0.25 


63 


2 


^'^ Kepler25c 


244.02 


6.24 


c 


349 


51 


3.98 


0.43 


127 


2 


^'^'^KeDler25b 


248.01 


7.20 


c 


368 


43 


8.87 


1 


108 


4 


^■2'^Kepler49b 


248.02 


10.91 


c 


354 


47 


16.9 


1.9 


66 


4 


^•^ Kcplcr49c 


250.01 


12.28 


c 


690 


260 


10.1 


1.1 


69 


4 


8 Koploi-26b 


250.02 


17.25 


c 


650 




8.1 


1.7 


45 


4 


^■^ Kepler26c 


262.01 


7.81 


c 


740 




25.2 


3.1 


108 


2 


5 Kepler50b 


262.02 


9.38 


c 


860 




16.1 


2.5 


87 


2 


^ Kepler50c 


270.01 


12.58 


c 


193.2 


10 


16.7 


5.2 


64 


2 


1 


271.02 


29.39 


c 


820 




13.2 


2.4 


22 


2 




277.01 


16.23 


c 


448 


21 


119.7 


3.3 


55 


1 


i'2'9Kepler36c 
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Table 4: - continued 



KOI 


Period 


Model 


Period 


ap 


Amp 




N 


Multi- 


Ref. 




[d] 




[d] 


[d] 


[min] 


[min] 




plicity 




308.01 


35.60 


C 


640 


79 


33.3 


2.3 


27 


1 


2 


314.02 


23.09 


C 


1200 




24.5 


2.5 


36 


3 




315.01 


35.59 


P 










25 


1 




318.01 


38.58 


c 


590 




4.57 


0.85 


19 


1 




319.01 


46.15 


c 


305 


18 


11.7 


1.4 


20 


1 




*341.01 


14.34 


c 


182 


11 


16.1 


2.7 


40 


2 




345.01 


29.88 


c 


820 




11.2 


1.9 


22 


1 




350.01 


12.99 


p 










55 


1 




377.01 


19.27 


c 


980 


200 


309.9 


3.2 


38 


3 


Kepler9b 


377.02 


38.85 


c 


1010 


180 


738.4 


6.5 


18 


3 


Kepler9c 


392.01 


33.42 


c 


1200 




67 


33 


23 


2 




410.01 


7.22 


p 










129 


1 




448.02 


43.59 


p 










19 


2 


2 


456.01 


13.70 


c 


720 




18.8 


1.7 


60 


2 


2 


473.01 


12.71 


c 


820 




30.9 


2.3 


60 


1 


2 


524.01 


4.59 


c 


331 


37 


17.5 


1.1 


196 


1 


2 


525.01 


11.53 


p 










83 


1 




*609.01 


4.40 


p 










208 


1 


16 


638.01 


23.64 


p 










23 


2 




738.01 


10.34 


c 


1600 




70 


11 


52 


2 


Kepler29b 


759.01 


32.63 


c 


1600 




73 


9.6 


22 


1 




760.01 


4.96 


p 










190 


1 




775.01 


16.38 


c 


212 


21 


15.9 


3.5 


40 


3 


5 Kepler52c 


775.02 


7.88 


c 


208 


14 


14.6 


2 


75 


3 


•"^ Kepler52b 
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Table 4: - continued 



KOI 


Period 


Model 


Period 


ap 


Amp 




N 


Multi- 


Ref. 




[d] 




[d] 


[d] 


[min] 


[min] 




plicity 




775.03 


36.45 


C 


305 


25 


11.6 


3.6 


15 


3 




784.01 


19.27 


C 


550 


110 


14.8 


4 


40 


2 


2 


806.02 


60.32 


c 


730 


290 


21.2 


1.6 


12 


3 


ii'i*epler30c 


829.01 


18.65 


c 


590 


150 


16.8 


3.5 


48 


3 


5 Kepler53b 


829.03 


38.56 


c 


474 


74 


25 


4.7 


23 


3 


^ Kepler53c 


841.01 


15.34 


c 


760 




14.2 


2.3 


51 


2 


Kepler27b 


854.01 


56.06 


c 


391 


42 


9.7 


4.1 


16 


1 


2 


872.01 


33.60 


c 


193 


6.9 


48.3 


5.1 


22 


2 


Kepler46b 


880.01 


26.44 


c 


940 




21.4 


3.6 


29 


4 




884.02 


20.49 


c 


805 


55 


178.6 


3.5 


36 


3 


1,5 


886.01 


8.01 


c 


1200 




99.6 


6.5 


78 


3 


Kcplcr54b 


886.02 


12.07 


c 


740 




23.2 


6.9 


36 


3 


^ Kepler54c 


902.01 


83.92 


p 










13 


1 




904.01 


2.21 


c 


256 


34 


4 


1.5 


196 


5 




904.02 


27.95 


p 










16 


5 


^ Kepler55b 


904.03 


42.15 


c 


800 




73 


11 


10 


5 


^ Kepler55c 


904.04 


4.62 


c 


196 


12 


7.6 


2.6 


97 


5 




904.05 


10.20 


c 


116.9 


3.7 


6.7 


2.8 


58 


5 




918.01 


39.64 


c 


1300 




14.6 


1.6 


23 


1 


2 


*935.01 


20.86 


c 


880 




19.5 


1.5 


38 


4 


Kepler31b 


*984.01 


4.29 


c 


487 


36 


46.4 


1.2 


62 


1 


2 


989.03 


16.16 


c 


1000 




47.8 


8.4 


36 


1 




1081.01 


9.96 


c 


1100 




89.9 


7.8 


85 


1 


2 


1102.01 


12.33 


c 


421 


50 


36.7 


5 


69 


2 


^'■^ Kepler24c 
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Table 4: - continued 



KOI 


Period 


Model 


Period 


crp 


Amp 




N 


Multi- 


Ref. 




[d] 




[d] 


[d] 


[min] 


[min] 




plicity 






1102.02 


8.15 


C 


416 


70 


24 


5 


93 


2 


2,7 


Kepler24b 


1145.01 


30.59 


C 


1600 




304 


14 


30 


1 


2 




1236.01 


35.74 


C 


760 


130 


68 


11 


20 


2 






*1241.01 


21.41 












37 


2 


5 


Kepler56c 


1241.02 


10.50 


C 


547 


100 


124 


24 


56 


2 


1,5 


Kepler56b 


1270.02 


11.61 


c 


455 


51 


26.3 


3 


43 


2 


2,5 


Kepler57c 


*1285.01 


0.94 












778 


1 


2,15 




*1452.01 


1.15 












767 


1 


15 




1529.01 


17.98 


c 


470 


110 


31.2 


7.6 


26 


2 


2,5 


Kepler59c 


*1546.01 


0.92 












675 


1 


15 




1573.01 


24.81 


c 


1500 




77.2 


3.8 


28 


1 


2 




1581.01 


29.54 


c 


730 




58 


14 


29 


1 


2 




1599.01 


20.42 


c 


770 




94 


19 


19 


1 


2 




1802.01 


5.25 


c 


240 


20 


7.85 


0.64 


183 


1 






1840.01 


7.04 


c 


1500 




67.4 


5 


125 


1 


2 




LS.lG.Ol 


4(i.:-')() 


c 


7!)0 




.S!)..") 


4 


21 
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Table 4: - continued 



KOI 


Period 


Model 


Period 


ap 


Amp 




N 


Multi- 


Ref. 




[d] 




[d] 


[d] 


[min] 


[min] 




plicity 




1884.01 


23.11 


P 










21 


2 




1907.01 


11.35 


C 


620 




10.1 


1.6 


72 


1 




2038.01 


8.31 


C 


980 




37.3 


5.1 


78 


4 





Ipll 

2F12 



^ Cochran et al 



^ Lissauer et al. 



Ford et al 



^ StefFen et al 



Ballard et al 



mm 



StefFen et al. 



mm 



mm 



( 2011a 



mm 



Carter et al 



( 2012b ) 



Holman et al 



mm 



mm 



^ Wabrvckv et al, 
^^ Nesvornv et al 



mm 



mm 



Wang et al 



mm 



Hxinglev et al 



Szabo et al 



mm 



(l2012bh 



Santerne et al 



mm 



* Discussed in Section [5] 
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Fig. 2. — KOIs with significant TTVs. The fit, plotted with a dashed hne, is composed of a 
straight hne together with a cosine function with the best found period. In one case — KOI- 
142.01, we fitted a straight hne with tvjo different cosine functions. When the periodicity 
could not be constrained, like in KOI-448.02, we fitted the data with a parabola only. In 
such cases, a note is written in the plot. 
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Fig. 3. — KOIs with significant TTVs. See Figure [2] for details. 
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Fig. 4. — KOIs with significant TTVs. See Figure [2] for details. 
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Fig. 5. — KOIs with significant TTVs. See Figure [2] for details. 
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Fig. 6. — KOIs with significant TTVs. See Figure [2] for details. 
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Fig. 7. — The KOIs with significant TTVs. See Figure |2] for details. 
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Fig. 8. — The KOIs with significant TTVs. See Figure |2] for details. 
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Fig. 9. — The KOIs with significant TTVs. See Figure |2] for details. 
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Fig. 10. — The KOIs with significant TTVs. See Figure |2] for details. 
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Fig. 11. — The KOIs with significant TTVs. See Figure |2]for details. 
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Fig. 12. — The KOIs with significant TTVs. See Figure |2] for details. 
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Fig. 13.— 



Time (BJD-2454900) 

The KOIs with significant TTVs. See Figure |2] for details. 
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5. KOIs with short-period TTV 

Our LS analysis of the transiting KOIs yields also 9 systems with highly significant 
short-period TTV in the range of 3 to 45 days. The modulation amplitudes are relatively 
small, in the range of 0.3 to 22 min, and it can be detected only because of the periodicity 
of the modulations and the long time span of the data. Figures [T1HT6] show the LS 
periodograms and the folded TTVs of the 9 systems, where one can see the prominent peaks 
of the periodograms. Table O lists the periods and amplitudes found. The table includes 
references to Section [6l where we briefly comment on these systems. 



As pointed out by ,Szab6 et al.l (|2012b|, hereafter Szl2), not all detected short-period 
modulations are due to physical TTVs. Apparent TTV periodicity can be induced either by 
the long cadence sampling of Kepler, or by interference with stellar periodic activity. 

The finite sampling rate of the observations can cause a shift of the orbital phases of 
the observations during a transit, inducing an apparent shift of the derived timing of that 
transit. This can evoke periodic TTV, with a period of 



Pi 



orb 



induced 



^ samp 



(3) 



where Porb is the orbital period, Psamp is the sampling cadence and [a;J is the floor of x 
(Szl2). Note that the induced periodicity is observed with "sampling" intervals equal to 
the planetary orbit, and therefore the Nyquist frequency of the observations is l/2Porb- If 
the frequency of Pmduced is larger, we would detect one of its aliases. We found two cases 
which suggest that the TTV included such effect (see below). 



The other effect is due to the stellar spot activity which modulates the stella r intensity 



with the stellar rotational period. Spot crossing (e.g.. 



Sanchis-Ojeda et al 



2OI2I) during a 



transit, or a different slope of the stellar brightness during a transit, can cause a shift in the 
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derived transit timing, inducing an apparent TTV periodicity with the stellar rotational 
period. In fact, all 9 systems show high level of stellar activity, and we had to check whether 
the detected TTV periodicity is due to that activity. 

To find the frequency of the presumed sampling-induced periodicity, we used for each 
KOI its Porb from Table [H and the pertinent Psamp, which was ~ 29.424 min for the long 



cadence, taken to be the median of the differe nces of the observed timings o: 



McQuillan. Aigrain fc Mazeh 



that KOI. We 



20121 ). and, if 



also searched for stellar spot periodicity (e.g. 
presented, checked if its frequency, or one of its aliases, was equal to the TTV one. We 
plotted the pertinent frequencies in Figures [HHTB1 



Table 5: KOIs with significant short period TTVs 



KOI 


Period 


Period 


ap 


Amp 


a-A 


N 


Multi- 


Ref. 




[d] 


[d] 


[d] 


[min] 


[min] 




plicity 




*258.01 


4.16 


36.01 


0.45 


7 


1.1 


175 


1 




*312.01 


11.58 


39.9 


0.48 


22.7 


1.9 


67 


2 




*725.01 


7.30 


22.07 


0.12 


5.67 


0.34 


129 


1 




*882.01 


1.96 


45.05 


0.61 


0.326 


0.075 


466 


1 


1 


*883.01 


2.69 


9.06 


0.045 


0.451 


0.046 


348 


1 


1 


*895.01 


4.41 


32.72 


0.37 


0.74 


0.12 


203 


1 


1 


*1128.01 


0.97 


3.36 


0.0023 


2.91 


0.52 


544 


1 




*1152.01 


4.72 


11.88 


0.042 


0.533 


0.06 


155 


1 


1 


*1382.01 


4.20 


34.38 


0.6 


1.109 


0.081 


194 


1 


1 



iSzl2 

* Discussed in Section |6] 
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Fig. 14. — The KOIs with short-period TTVs. For each KOI, the LS periodogram and 
the folded TTV are plotted. The dashed green line represents the stellar activity frequency 
or one of its aliases, if presented in the stellar light curve, and the dash-dotted red line 
represents the frequency induced by the sampling. 




Fig. 15. — The KOIs with short-period TTVs. For each KOI, the LS periodogram and the 
folded TTV are plotted (see Figure [TH for details). 
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Fig. 16. — The KOIs with short-period TTVs. For each KOI, the LS periodogram and the 
folded TTV are plotted (see Figure [TH for details). 
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6. Comments on Individual Systems 

In this section we comment on a few KOIs from Tables HHS] In particular, we folded 
the light curves of all 100 systems with their orbital period, and looked for a secondary 
eclipse. For eleven systems we find significant secondary dip, in most cases at phase ~ 0.5, 
which we interpret as either an eclipse of a secondary star, or as a planetary occultation. 
In one case we believe that the published period is twice the correct one. We also point 
out when we suspect that the periodic TTV could be induced either by the long cadence 
sampling or by the stellar spot periodic activity. Some of these systems were analyzed in a 
similar way by Szl2, who used six quarters of Kepler data to look for TTV periodicity. 



KOI-142.01 (Figure 2): The TTV modulation, with one of the largest amplitudes, 
could not be fitted with one cosine function. The TTV includes at least two different 
frequencies, as can be seen in its figure. This could be induced by two additional 
planets. 



KOI-190.01 (Figure 3): ISanterne et al.l ( 120121 ) found this system to be a diluted 
eclipsing binary. 

KOI-258.01 (Figure 14): The folded light curve reveals a dip at phase ~ 0.5, probably 
a secondary eclipse, and therefore the system is probably an eclipsing binary (=EB). 
The star has significantly high level of activity, probably due to stellar pulsations. 
The TTV LS periodogram displays two prominent peaks, one is the first harmonic of 
the other. 

KOI-312.01 (Figure 14): The autocorrelation of the stellar photometry reveals a weak 
but stable modulation with a short period of 0.17073 d. The green line in the figure 
is an alias of this frequency. 
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KOI-341.01 (Figure 5): The orbital period is probably half the one published and 
used here. 



KOI-609.01 (Figure 7): ISanterne et al.l ( 120121 ) found this system to be a diluted 



eclipsing binary. The folded light curve reveals a dip at phase ~ 0.5. 

KOI-725.01 (Figure 14): The folded light curve displays a shallow occultation. The 
stellar photometry shows a pulsation of 8.58 d. The TTV LS periodogram displays 
two prominent peaks, one is the first harmonic of the other. 

KOI-882.01 (Figure 15): The photometry displays strong stellar pulsations with a 
frequency of 3.921 d, very close to twice the orbital period. The second peak in the 
periodogram is an alias of the first one, relative to the pulsation Nyquist frequency. 
Szl2 found the same TTV periodicity, with noisier periodogram. 

KOI-883.01 (Figure 15): The TTV LS periodogram displays two prominent peaks. 
The higher one coincides with the stellar rotation, with a period of 9.02 d, and the 
smaller one with the induced sampling frequency. Szl2 reached the same conclusion. 

KOI-895.01 (Figure 15): The folded light curve displays a shallow occultation. The 
star has very pronounced activity, with a period of 5.0935 d, probably due to stellar 
spots. The TTV LS periodogram displays one prominent peak, with frequency that 
coincides with the beat frequency with the stellar one. Szl2 identified a different TTV 
period, close to the third harmonic of the frequency we find. 

KOI-935.01 (Figure 10): The folded light curve probably displays an occultation. 

KOI-984.01 (Figure 10): TTV deviates from the strictly cosine function at 
BJD ~ 2454900 + 100. 
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• KOI-1 128.01 (Figure 16): The TTV LS periodogram displays one prominent peak, 
with frequency that coincides with the induced samphng one. 

• KOI-1 152.01 (Figure 16): The TTV LS periodogram displays one prominent peak, 
with frequency very close to the stellar rotational one, at a period of 2.95 d. Szl2 
found a TTV period which is twice the period we found, and did not associated it 
with the stellar modulation. They also detected a secondary echpse not in phase 0.5, 
and concluded that the system is eccentric EB. 

• KOI-1241.01 (Figure 11): This planet is part of a known multiple system. KOI-1241.02 

shows a clear periodic TTV, and KOl-1241.01 might have anti correlated TTV, as 
expected. However, the significance of the periodicity of 1241.01 is not high enough 
to be included in Table 4. 

• KOI-1285.01 (Figure 11): The folded light curve, with the orbital period of 0.9374 d, 
reveals a dip at phase ~ 0.5, probably a secondary eclipse, and therefore the system 
is probably an EB. The TTVs display coherent modulations, but not a clear stable 
periodicity. The star has a significant high level of periodic activity, with a period of 
0.9362 d, which is close but not identical to the orbital period. Szl2 identified a few 
different possible TTV periods. They suspected that two of their periods are affected 
by the stellar modulation. 

• KOI-1382.01 (Figure 16): The folded hght curve reveals a dip at phase ~ 0.5, 
probably a secondary eclipse, and therefore the system is probably an EB. The TTV 
LS periodogram displays one prominent peak, with a frequency that coincides with 
one of the ahases of the stellar rotational one, at a period of 4.79 d. Szl2 identified 
the same TTV period, although with much higher first harmonic. They failed to 
notice that the TTV periodicity comes from the stellar rotation. 
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• KOI-1452.01 (Figure 11): The folded light curve, with the orbital period of 1.1522 d, 
reveals a dip at phase ~ 0.5, probably a secondary eclipse, and therefore the system 
is probably an EB. The TTVs display coherent modulations, but not a clear stable 
periodicity. The stellar photometry displays strong stellar pulsations with frequencies 
of 0.65597, 0.7097 and 0.83 d^^ Szl2 found it to be a multi-periodic candidate. 

• KOI-1546.01 (Figure 12): The folded hght curve, with the orbital period of 0.9176 d, 
reveals a dip at phase ~ 0.5, probably a secondary eclipse, and therefore the system is 
probably an EB. The TTVs display coherent modulations until BJD ~ 2454900 + 600, 
but not a clear stable periodicity. The stellar photometry displays a strong periodicity 
of 0.933 d, probably due to stellar rotation, a period very close but not identical to 
the orbital period. Szl2 found it to be a multi-periodic candidate. 



7. Discussion 

We presented here 100 KOIs with highly significant TTVs, 91 with long-term 
modulations (Section 4, Table 4), and nine KOIs with short-period low-amplitude TTV 
periodicities (Section 5, Table 5). Out of the 91 systems, 76 showed clear periodicities, 
with well determined periods and amplitudes. Another 11 KOIs had periods too long to 
be established without a doubt. For those we need to wait for more data before the TTV 
period can be safely determined. Another four systems display coherent modulations, but 
not a clear stable periodicity. 

As discussed in Section 6, some TTVs, the nine short-period ones in particular, can be 
due to some interference with stellar photometric modulation, or some other "instrumental" 
effect. However, let us assume that most of the observed long-term periodic TTVs of 
the Kepler KOIs are caused by dynamical interaction between planets, as predicted by 
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Agol et al.l ( 120051 ) and lHolman fc Murrayl (|2005[ ). For the single-planet KOIs, let us further 
assume that the interaction is induced by an unseen additional planet. We did not detect 
the transits of the additional planet either because of a larger inclination angle relative to 
our line of sight, or because its radius was too small. We presented here 37 single-KOI 
planets with significant TTVs, 7 of which were with TTV periodicity probably longer than 
~ 2000 days. On the other hand, only 4 out of 54 multiple KOIs with significant TTVs 
showed an indication for long-term periodicity. 

The different fraction of KOIs with long-term modulations among the single and 
multiple KOIs, if real, is intriguing. The longer TTV period could be due to the fact that 
the additional unseen planet is closer to a mean-motion resonance. If this is true, we can 
conjecture that pairs of planets with orbital periods close to a resonance tend to have larger 
relative inclinations. Another possibility is that we did not detect the additional planets in 
the single-KOI systems because the other planet is either too small or has a large orbit. 
Both features are expected to produce longer TTV periodicities. 

We also find an indication for some correlation between the KOI period and the period 
of its TTV, as can be seen in Figure [T71 This is of no surprise, as the orbital period of a 
planet is the natural time unit of the dynamical interaction, and therefore we expect the 
TTV periodicity to be correlated with this time unit. Another correlation, between the 
amplitudes and the periods of the detected TTV periodici ties, emerges from the present 
sample (see Figure [T8|) . as was predicted, for example, by 



Agol et al 



(|2005|). The same 



correlation appears when we plot the amplitude in units of the KOI orbital period. 

We pointed out possible non-physical origin of some of the TTVs presented here. 
In particular, the short-period modulations could be due to either to the long cadence 
sampling of Kepler oi by stellar spot periodic activity (Szl2). We found evidence that 4 out 
of the 9 short-period detected TTVs are due to the stellar periodicity. We also found that 
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TTV period as a function of orbital period 



A multiples 
• singles 



10 



T3 
O 



CD 
Q. 



% 

A 
A 



^ A • 

A • A 



A 
A 



Orbital period [d] 

Fig. 17. — The TTV period as a function of the orbital period. Plot is on a log- log scale. 
Single KOIs are plotted with red circles, and multiples with green triangles. 
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Sine amplitude as a function of the TTV period 
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Fig. 18. — The amplitude of the TTV modulation as a function of its period. Plot is on a 
log-log scale. Single KOIs are plotted with red circles, and multiples with green triangles. 



-44- 



KOI-883.01 and 1128.01 show a periodicity induced by the Kepler samphng. 



The sample of 100 KOIs with significant TTVs includes 54 systems discussed by Fll 
(16 KOIs), F12 (36 KOIs) and S12 (18 KOIs), all based only on a fraction of the data 
available now. References to those three works can be found in Table |H It is interesting 
to compare the analysis of Fll, F12 and S12 on one hand and the present results on the 
other hand, and see how doubling the time span can change our assessment of the nature 
of the modulation. In many cases the time span of the first 6 quarters was not long enough 
to detect a local maximum and minimum of the TTV modulation, and therefore the 
periodicity of the modulation could not be estimated. One illustrious example is KOI-142, 
with its peak-to-peak amplitude of more than 1200 min. 

One could hope that the accumulating details of the observed TTV could give some 
hints for the orbital elements of the pert urbing unseen plane t , at le ast f or some of the sing le 



KOIs. However, as discussed already by 



Holman fc Murray! (120051 ) and 



Agol et al 



(120051) 



the amplitude and periodicity of the TTV modulation depends on various parameters, in 
particular the mass and orbital period of the unseen planet and how close the orbits of the 
two planets are to some mean motion resonance. Therefore, it is quite difficult to deduce 
the parameters of the unseen planet, although some stringent constraints can be derived. 



as was done by 



Ballard et al. 



(I2OIII ) for Kepler-19. We hope that the available catalog will 



push towards similar work on other single-KOI systems with significant TTV. 

One parameter that bears interesting implications on our understanding of the 
planetary formation is the relative inclination between the plane of motion of the observed 
planet and that of the presumed interacting planet for the single KOI with significant TTV. 
Relative inclination can induce a precession of the orbital motion of the observed planet, 
which can manifest itself in a modulation of the transit duration and depth (TDV and 
TPV). Although the focus of the present work is on the TTVs, the catalog, which includes 
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derived TDV and TPV, can, in principle, help to point out systems with relative inclination. 
Furthermore, stringent upper limit on TDV and TPV for systems with detected TTVs can 
help to constrain the relative inclination between the planets. 

However, an observed precession is not necessarily induced by a planet with a 
non- vanishing relative inclination. An observed precession just proves that the total angular 
momentum of the system is not parallel to the orbital angular momentum of the transiting 
planet. The origin of the precession could be misalignment of the stellar rotation axis 
relative to the angular momentum of the planet, an idea t hat was unt hinkable not long ago, 
but has now sol id evidence in the accumulating data (see 



Winn 



idea is KOI-13 flSzabo et al. 



2011 



201 ll ). An example for this 



see also Figured]). Regardless, we suggest that systems 



with detected significant TDV and TPV deserve further close study. 

Finally, we present here a systematic TTV analysis of twelve quarters of Kepler 
observations of all KOIs. One could expect that the derived TTVs, for the single KOIs in 
particular, could help constructing a statistical picture of the frequency an d architecture o: 



the p opulation of the planetary multiple systems of the Kepler KOIs (e.g.. 



Ford et al 



2011 



20121 ) , like the discussion that followed the disc overies of if ep/er multiple transiting planets 



le.g- 



Lissauer et al. 



2011b 



Steffen et al. 



20101 ). To perform such a statistical analysis 



one needs to model the dependence of the detectability of TTV coherent modulation on 
the parameters of the unseen perturbing planet, and to convolve it with the observed 
parameters of the transit, like its SNR and duration, together with the orbital period of the 
transiting planet. The present catalog can be used for such a study. 



The research leading to these results has received funding from the European Research 
Council under the EU's Seventh Framework Programme (FP7/(2007-2013)/ ERG Grant 
Agreement No. 291352). All photometric data presented in this paper were obtained from 
the Mikulsky Archive for Space Telescopes (MAST). STScI is operated by the Association of 



-46 - 



Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. Support 
for MAST for non-HST data is provided by the NASA Office of Space Science via grant 
NNX09AF08G and by other grants and contracts. 



-47- 

REFERENCES 

Agol, E., Steffen, J., Sari, R., & Clarkson, W. 2005, MNRAS, 359, 567 
Batalha, N. M., Rowe, J. F., Bryson, S. T., et al. 2012. larXiv: 1202.58521 
Ballard, S., Fabrycky, D., Fressin, F., et al. 2011, ApJ, 743, 200 
Borucki, W. J., Koch, D. G., Basri, G., et al. 2011, ApJ, 736, 19 
Carter, J. A., Agol, E., Chaplin, W. J., et al. 2012, Science, 337, 556 
Cochran, W. D., Fabrycky, D. C, Torres, G., et al. 2011, ApJS, 197, 7 
Fabrycky, D. C, Ford, E. B., Steffen, J. H., et al. 2012, ApJ, 750, 114 
Ford, E. B., Rowe, J. F., Fabrycky, D. C, et al. 2011, ApJS, 197, 2 
Ford, E. B., Fabrycky, D. C, Steffen, J. H., et al. 2012, ApJ, 750, 113 
Ford, E. B., Ragozzine, D., Rowe, J. F., et al. 2012, ApJ, 756, 185 
Holman, M. J., & Murray, N. W. 2005, Science, 307, 1288 
Holman, M. J., Fabrycky, D. C, Ragozzine, D., et al. 2010, Science, 330, 51 
Lissauer, J. J., Fabrycky, D. C, Ford, E. B., et al. 2011a, Nature, 470, 53 
Lissauer, J. J., Ragozzine, D., Fabrycky, D. C, et al. 2011b, ApJS, 197, 8 
McQuillan, A., Aigrain, S., & Mazeh, T. 2012, submitted to MNRAS 
Mandel, K., & Agol, E. 2002, ApJ, 580, L171 

Mazeh, T., Nachmani, G., Sokol, G., Faigler, S., & Zucker, S. 2012, A&A, 541, A56 
Nesvorny, D., Kipping, D. M., Buchhave, L. A., et al. 2012, Science, 336, 1133 



-48 - 

Sanchis-Ojeda, R., Fabrycky, D. C, Winn, J. N., et al. 2012, Nature, 487, 449 

Santerne, A., Diaz, R. F., Moutou, C, et al. 2012, A&A, 545, A76 

Steffen, J. H., Batalha, N. M., Borucki, W. J., et al. 2010, ApJ, 725, 1226 

Steffen, J. H., Ford, E. B., Rowe, J. F., et al. 2012, ApJ, 756, 186 

Steffen, J. H., Fabrycky, D. C, Ford, E. B., et al. 2012, MNRAS, 421, 2342 

Steffen, J. H., Fabrycky, D. C, Agol, E., et al. 2013, MNRAS, 428, 1077 

Szabo, G. M., Szabo, R., Benko, J. M., et al. 2011, ApJ, 736, L4 

Szabo, G. M., Pal, A., Derekas, A., et al. 2012, MNRAS, 421, L122 

Szabo, R., Szabo, G. M., & Dalya, G. 2012. larXiv:1207. 72291 

Tamuz, O., Mazeh, T., & North, R 2006, MNRAS, 367, 1521 

Tingley, B., Palle, E., Parviainen, H., et al. 2011, A&A, 536, L9 

Wang, S., Ji, J., & Zhou, J.-L. 2012, ApJ, 753, 170 

Winn, J. N. 2011, European Physical Journal Web of Conferences, 11, 5002 



This manuscript was prepared with the A AS lATgX macros v5.2. 



